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Abstract. An extensive search for the complex molecule trans-ethyl methyl ether towards several hot core regions 
has been performed. Using the IRAM 30m telescope and the SEST 15m we looked at several frequencies where 
trans-ethyl methyl ether has strong transitions, as well as lines which are particularly sensitive to the physical 
conditions in which the molecule can be found. We included G34.26, NGC6334(I), Orion KL, and W51e2 which 
have previously been proven to have a rich chemistry of complex molecules. Our observations cannot confirm the 
tentative Orion KL detection made by |Charnley et al.(2001)] within their stated column density limits, but we 
confirm the existence of the trans-ethyl methyl ether towards W51e2 with a column density of 2 x 10 14 cm -2 . 
The dimethyl ether/methanol ratio of 0.6 as well as the newly found ethyl methyl ether/ethanol ratio of 0.13 
indicate relative high abundances of ethers toward W51e2. Furthermore, the observation of ethyl methyl ether 
also confirms the importance of ethanol as a grain mantle constituent. We present new upper limits of around 
8xl0 13 cm -2 for the column densities of the molecule toward Orion KL, G34.26, NGC6334(I) and estimate the 
column density towards SgrB2(N) to be of the same order. The W51e2 observations are discussed in more detail. 
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1. Introduction 

Various large and complex molecules of prebiotic impor- 
tance have been found as constituents of dense interstel- 
lar clouds but their detection remains difficult and some 
of them doubtful Snyder et al.(2005) . Fortunately, space 



offers many places with a rich chemistry and hot molec- 
ular cores (HMCs) are prime targets for such investiga- 
tions. However, the spectra towards hot cores often reveal 
a high line density and thus proper line assignments are 
difficult and have to be viewed with a critical eye. The 
ethyl methyl ether (CH3OC2H5) is an asymmetric top 
molecule with two methyl groups acting as internal rotors. 
Ethyl methyl ether (EME) has two conformers where the 
trans-EME is lower in energy than the gauche-EME by 
550 cm -1 Durig et al.(2002')l and thus favored as an as- 
trophysically detectable species. We searched for several 
transitions of trans-EME to confirm the tentative detec- 



tion of this molecule by Charnley et al.(2001) A positive 
result means the detection of one of the biggest molecules 
in the interstellar medium (ISM). Charnley et al.(2001) 
observed a single line in W51 el/e2 and Orion KL at 
160.1 GHz and another line in SgrB2(N) at 79.6 GHz 
which they assigned to trans-ethyl methyl ether. Based 
on these isolated measurements, Charnley et al.(2001) es- 



timate the column density of EME to be in the range of 
10 14 - 10 15 cm -2 in Sgr B2(N). Thus, assuming a hydrogen 
column density of 10 24 cm -2 gives a fractional abundance 
of lO" 10 - 10" 9 . For W51, they estimate EME to have 
an abundance of a few times 10 -10 . Furthermore, they 
conclude that the emission of EME closely follows that 
of methanol, as expected if these molecules are chemi- 
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cally linked. Charnley et al.(1995) suggest the following 
production mechanism for EME in hot cores where tem- 
peratures can rise to 100-300 K. Icy mantles which are 
thought to contain alcohols (e.g. methanol, ethanol) evap- 
orate from dusty grains, e.g. Millar & Hatchell(1998) the 
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alcohols themselves are thought to be produced via surface 
reactions starting from the precursor CO. Once the man- 
tles evaporate the alcohols can be converted into ethers 
and other species. The production of ethers is thought to 
happen via ion- molecule reactions 1 (i), methylisation of 
the protonated species 2 (ii) and dissociative recombina- 
tion with electrons (iii): 



i) C 2 H 5 OH C 2 H 5 OH 



HCO+ 



ii) C 2 H 5 OH^ 



CH3OC2H+ 



CH3OH 

iii) CH 3 OC 2 H+ — CH 3 OC 2 H 5 

In hot cores saturated species are the dominant con- 
stituents. For instance, methanol CH3OH is highly 
abundant and dimethyl ether (CH^O and ethanol 
C2H5OH are seen only in hot cores |Herbst(2001)"l . 
Charnley et al.(1995) mention in their work that dimethyl 
ether, diethyl ether (C 2 H 5 ) 2 0, and trans-EME should 
be present in detectable quantities within cores rich in 
ethanol and methanol. Furthermore, a detection of either 
EME or diethyl ether would confirm the importance of 
ethanol as a grain mantle constituent. 



2. EME Spectrum and Observations 

Only two investigations on the laboratory rotational spec- 
trum of EME have been reported in the literature. The 



more recent measurements by Fuchs et al.(2003) from 55 
GHz to 350 GHz also include the previously published 



data below 40 GHz by Hayashi fc Kuwada(1975)| in their 
analysis. Reliable frequency predictions up to 400 GHz are 
now available. Since EME consists of two methyl groups 
designated as internal rotors, the laboratory spectra re- 
veal extra line splittings which are not necessarily seen 
in space due to larger line widths in the HMC environ- 
ment. The spectroscopic aspects of ethyl methyl ether in 
the mm-wave region are described in our earlier publica- 
tions (see Fuchs et al.(2003)| Fuchs(2003) and references 
therein). EME is an asymmetric top molecule which has 
weak a-type transitions and fairly strong &-type transi- 
tions above 100 GHz which in theory split into quintets 
caused by internal torsion. Laboratory measurements can 
resolve only doublets and quartets with splittings mostly 
between 0.8 and 2 MHz in our examined frequency ranges. 
However, none of the astronomically observed lines is tor- 
sionally resolved. Fig.^shows the complexity of the EME 
spectrum up to 400 GHz with its characteristic features. 

The IRAM 30m telescope at Pico Veleta, Spain, has 
been used during several days in Sept 2002 and July 2003 
and the SEST 15m telescope at La Silla, Chile, in March 
2003. The hot core regions of G34.26, NGC6334I, Orion 
KL and W51e2 were observed, see Table Rotational 
transitions of EME have been searched in the 1mm-, 2mm- 
, and 3mm-range. Depending on the source, between 5 



see 
see 



Mackay et al.(1979) 



Karpas & Meot-Ner 



|Bohme et al. 
(1989)| 



(1979) 



and 11 frequency settings have been used, see Table [S] for 
observational details. Typical frequency resolutions were 
0.3, 1 and 1.25 MHz, typical bandwidths between 490 and 
1020 MHz and rms noise levels between 0.01 K at 81 GHz 
and 0.02 - 0.05 K at 245 GHz. The bands were chosen 
to allow for the observation of CH3OH (a-type), C2H5OH 
and CH3OCH3 transitions as well. At least 5 strong lines 
have been measured from each of these species, except for 
CH3OH in Orion KL where only three transitions have 
been observed. 

The intensities are measured in T£ the effective an- 
tenna temperature corrected for spillover losses and at- 
mospheric attenuation, and are then converted to T m b = 
T^/p using the beam efficiencies p of the telescopes at the 
given frequencies. 

All sources were examined using long integration times 
up to 220 minutes on+off time, elevations greater 40° and 
under good atmospheric weather conditions. 

3. Source selection 

A selection of previous work on the sources is given below. 
The authors do not intend to give a review of these sources 
but rather focus on contributions relevant for the EME 
analysis. 

G34.3+0.15 has be en surveyed by 
and 



MacDonald et al.(1996) 



in the 330-360 GHz 



Thompson et al.(1999) 



Nomura fc Millar(2004) 



region 

examined the chemical structure of this region and gave 
an overview of the observed species and their abundances 
along with some comparisons to their model predictions. 
Single molecule investigations like ethanol have been 



performed by Millar et al.(1995) Nummelin et al . (1998b) 
examined CH3CHO and c-C 2 H40 but also gave abun- 
dances for methanol, ethanol and dimethyl ether (DME). 
Hcaton et al.(1989) mapped the G34.3 region using the 
NHs(3,3) emission and determined a size of 3" x 1".2. 
The column density of H 2 towards G34.26 has been 
estimated b y |Millar et al.(1995)| to be N t (H 2 ) = 5.3 x 
10 23 cm -2 . Nummelin et al. (1998b) published a similar 
result of JVt(H 2 ) = 3.0 x 10^ 3 cm" 2 for a 20" source size. 

NGC6334I (NGC6334F, G351.4 1+0.64) has 
been 



mapped 



by Sandell(2000) 



McCutcheon et al.(2000) took several spectra in the 
range between 334 and 348 GHz using the JCMT. 
They also mapped the region in CH3OCH3 emission. 
Nummelin et al. (1998b) estimated the abundances of 
ethanol, methanol and DME in NGC6334I and a hydro- 
gen column density of JV t (H 2 ) = 2.0 x 10 23 cm~ 2 within 
20". |Cheung et al.(1978)| estimate N Ha » 10 24 cur 2 . 
Orion KL is a well studied object in the 



sky. |Turner(1989)[ Turner(1991) performed a line sur 



vey between 70 and 115 GHz 
sented line spectra between 159.7 to 164.7 GHz 



Lee fc Cho(2002) 



pre- 
and 

Sutton et al. (1985)| scanned the spectral region between 
215 and 247 GHz. |Blake et al.(1987)| studied the chem- 
ical composition of OMC-1 in the 230 GHz atmo- 
spheric window. They estimate the hydrogen column den- 
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sity to be N# 2 = 1.0 x 10 3 cm for the hot core. 
The region between 330 to 355 GHz was examined by 
Sutton et al.(199lJt and Schilke et al.(1997) performed a 
survey from 325 GHz to 360 GHz ( sec also references 
therein). The distribution of some key molecules for the 



EME analysis has been presented by Sutton et al.(1995) 



Caselli et al.(1993) applied a dynamical-chemical model 
which included abundances of many observed species. 
The W51 molecular cloud is described in 



Carpenter fc Sanders(1998)| The chemistry of the W51e2 
region has been investigated by several groups. For ex- 
ample, |Zhang et al.(1998)]has done CH 3 CN and CS (3-2) 



observations, Remijan et al.(2004) Remijan et al.(2002) 



investigated CH 3 CN and CH 3 COOT I~lLiu et al.(200"T) 



looked for HCOOH in this source and Kuan et al.(2004) 
searched for several molecules of potential prebiotic 
importance. The hydrogen column density 3.6 x 10 23 
cm -2 Ikeda et al.(2001) is based on C 18 data from 
Schloerb et al.(1987)| 
Also the SgrB2 
many groups. It is 
tions but since Charnley et al.(2001) detected one line 
at an EME transition frequency we use the pub- 
lished data to estimate the column density of EME 
in Section A spectral line survey from 70 and 150 



region has been 
not part of 



investigated by 
our own observa- 



GHz has been performed by Cummins et al.(1986)[ and 
from 70 to 115 GHz by |Turner(1989)| |Turner( 1991)1 
A further 3mm spectral line survey has been carried 
by |Friedel et al.(2004)| using the NRAQ12m and 



out 

the BIMA array telescopes. Nummelin et al. (1998a) ex- 
tended these observations to frequencies between 218 and 
263 GHz. Selected molecules, like the oxygen bearing 
CH3COOH, have been studied by |Remijan et al.(2002)| 
and NH 3 by Hiittemeister et al.(1993) with comparisons 
between several sources in cluding those studied here. 
Marti'n-Pintado et al.(1990)| estimated N ff2 = 5.0 x 10 24 
cm" 



2 to be the hydrogen column density 3 
Ikeda et al (2001)] 



|Ikeda et al.(2002~ gave temperatures and abundances 
of methanol, ethanol and in most cases also for DME for all 
above mentioned sources. Table|21summerizes the CH 3 OH, 
C2H5OH, CH3OCH3 and C2H5CH3 abundances including 
our latest results. 



4. Data Analysis method 

For the analysis of the observed lines we used the Grenoble 
Image and Line Data Analysis Software GILDAS 4 and in 
addition the myXCLASS code written by Peter Schilke 
which is backed by a database of astrophysically impor- 
tant molecules, including the EME data. myXCLASS can 
simulate the spectra of many molecules at a time, includ- 
ing their line profiles within the GILDAS environment 
and has restricted fitting capabilities. Fig. |2] shows the 



cm 

4 



Note, |Charnley et al.(2001)| assume Nh 2 = 1.0 x 10 
see |http: / / www.iram.fr/IRAMFR/GILDASI 



observed spectrum of W51e2 at 151 GHz and the syn- 
thetic spectrum calculated by myXCLASS assuming op- 
tically thin lines where for each molecule the excitation 
is characterized by a rotational thermal equilibrium tem- 
perature. For each source a simulated spectrum has been 
created which covered all measured frequency bands. This 
enables an estimate of the abundances and temperatures 
of the molecules in these sources and can help to reveal 
blendings of several lines at a given frequency position. 
myXCLASS is specially suited to deal with many convo- 
luting line profiles with different line widths and intensi- 
ties as can regularly be observed in hot core regions where 
lines of different molecules blend with each other. 

The spatial extent of the molecule can influence the 
value of the inferred column density. This is because the 
circular Gaussian beam size diameter 9 b of the telescope 
varies with frequency v, i.e. 9b ~ l/v (see telescope 
beam HPBW in Table EJ. If the telescope beam points 
towards a source of circular Gaussian type with size 9s 
the resulting convolution beam filling factor b is given by 
b = 9l/(9 2 B + 91), seelSnyder et al.(2005)l and the work 
by Ulich & Haas(1976) In the case of 9s S> 9b the beam 
filling factor is close to unity, for 9s ~ 9b the beam dilu- 
tion is 0.5 and for 9 S < 9 B it follows b « 9 S /9 2 B . Thus, 
assuming that the intrinsic intensity distribution Tg of a 
source is Gaussian, a beam filling factor can be intro- 
duced to calculate the measured intensity T m b = b ■ Tg. 
HMCs are compact objects with diameters < 10". Using 
the IRAM 30m telescope with HPBW around 30" and 
10" at 80 and 250 GHz respectively and assuming a 10" 
source size the beam filling factor is 6(@80GHz)= 0.1 and 
6(@250GHz)=0.5. Hence, an intensity difference of a fac- 
tor 5 can occur by not correcting for beam dilution. The 
program myXCLASS has the source size as an input pa- 
rameter and can thus simulate the effect of the assumed 
source size on the line intensities. 

For the EME lines we can expect small optical depths 
so that the rotational temperature T rot and the total 
molecular column density Ni [cm -2 ] can be determined 
by a least-squares fit of 



log 



/ 3fc B / T mb dv 



10 



log 



in 



Q(Trot) 



\og w {e)E u 

T-mt 



(i) 

with J T m i,dv the observed line integral, n 2 S the dipole 
and line strength value derived from Pickett et al.(1998) 
(see their Eq.(2)), Q the partition function, and E u 
the energy of the upper level. The resulting fit can be 
plotted in a rotational-temperature-diagram (RTD) or 
Boltzmann plot. For detections at the confusion limit 



Snyder et al.(2005) has shown that Boltzmann plots can- 
not necessarily be used as independent tools to verify 
interstellar spectral line identifications. Any false assign- 
ment does not automatically reveal itself in the RTD by 
large displacements from the fitted line and thus can only 
be avoided by careful line check in catalogs and the lit- 
erature. As a standard procedure any candidate line of 
EME has to be checked against competing lines, i.e. the 
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expected or known intensities of these other lines at the 
EME frequencies have to be determined. Therfore, a con- 
sistent and comprehensive list of molecule abundances in 
the source is important. Line surveys are invaluable aids 
in determining the chemical composition of a source. 

5. Results 

5.1. W51e2 

Several lines at EME transitions have been observed to- 
wards W51e2, see Fig. |3J We looked for strong EME fa- 
type lines and at the same time tried to exclude as many 
blendings with other molecules. By using myXCLASS we 
first estimated a column density of ~10 14 cm -2 and a ro- 
tational temperature of roughly 100K. We checked every 
frequency band for consistency by comparing a simulated 
spectrum of EME and other known molecules with the 
observed spectrum. Fig. 0] shows an observed line with 
simulated EME lines superimposed. The laboratory mea- 
surements of the EME line positions are at a precision level 
unmatched by any astronomical observations. The figure 
illustrates the effect of the linewidth on the total observ- 
able line shape and shows that the unique EME feature 
which is due to the torsional splitting remains unresolved 
in the observed spectrum. The observed lines 5 at transi- 
tion frequencies of EME are shown in Table |3J The more 
stronger lines can be found above 150 GHz whereas many 
lines at the la - 2a level have been observed at lower fre- 
quencies. Hot cores reveal a great richness and density in 
molecular emission lines, see Fig. |H1 and it is not aston- 
ishing that 77% of the EME lines in the observed bands 
arc blended by strong lines of other known molecules. In 
our case by far the most blendings occur in the 215 GHz 
band which is therefore not further detailed in Table 01 
12% of the blendings are at line positions of yet uniden- 
tified species. However, from intensity considerations it 
is clear that these latter lines cannot be dominated by 
EME because in this case EME should have also strong 
lines at other frequencies where we have not observed any 
strong feature. Note that many transition in TableOlhave 
labels with possible blendings. By carefully checking the 
databases and literature but also by fitting many other 
species to our spectra we could estimate the abundances 
of several species which could be potential competitors to 
the EME lines. To further clarify our detection we made 
two Boltzmann fits, as shown in Fig.|3 Here, every line in 
Tablc|31which is marked with an (•) or (□) is depicted. If 
we use all these 14 lines in the Boltzmann fit (I) we get 
T rot = 126±115 K and N t = 7.9±f x 10 14 cm~ 2 as the 
beam averaged column density. Since we have many weak 
lines (i.e. around la level) the reliability of the fit can be 
checked by excluding the weaker or doubtful lines marked 
with (□=noisy-like). Fit (II) only includes lines (•) with 
intensities larger than 2a. As a result we obtain T rot = 



54±38 K and N t = 6.3±f x 10 14 cm" 2 . This fit II should 
be regarded as more reliable than fit I. Furthermore, it has 
been shown that the excluded weak and noisy lines (□) 
are not in contradiction to fit II. Due to the large errors of 
the weak lines in the Boltzmann plot this is indeed not the 
case, but we would still expect slightly higher intensities 
at the (□) frequency positions. 



W51e2 is known to have a size < 10" so that the ef- 
fect of beam dilution should be considered. For each fit 
(I and II) the above mentioned temperatures and abun- 
dances of EME can be recalculated by assuming an 5" 
source size. This results in source averaged column den- 
sities and temperatures: fit I gives T rot = 74±50 K and 



1+29 



2Ii.5 x 10 cm ~ . fit 11 T iot = 69±94 K and N t = 
x 10 14 cm -2 . Hence, the column densities and tem- 
peratures of the diluted and non-diluted beam yield com- 
parable results within their error limits and we will use the 
source averaged value of fit (II) for further investigations. 
As has already been stated, RTDs cannot necessarily be 
used as reliable independent evidence to support the de- 
tection of a molecule close to the confusion limit, as was 



shown by Snyder et al.(2005) in the case of glycine. Our 
estimated column density is also consistent with obser- 
vations in less favored frequency regions as can be seen 
in Fig. HJ1 where an EME line is expected. For our fitted 
2 x 10 14 cm -2 column density no contradiction within the 
la level can be found, but for a N t = 4 - 8 x 10 14 cm~ 2 the 
simulated spectrum is in clear conflict with observation. 

It is commonly believed that most complex molecules 
exist within a small area around the forming star, 
e.g. around 5" for W51c2. On the other hand 



5 The 53,2-42,3 transition line at 160.1 GHz from 
|Charnley et al.(2001)| is not listed in Table Eland has not been 
considered in our fits. 



Hollis et al.(2001) have shown that complex molecules 
such as glycolaldehyde can be found on spatial scales of > 
60" in the case of SgrB2. With our data it is not possible 
to determine the size of the emitting region but the source 
size has surely an influence on the detected signal strength 
which we tried to estimate using the myXCLASS program. 
Table 0] lists the result of possible column densities con- 
sistent with the observations if we assume certain tem- 
peratures between 30 K and 150 K. |Charnley et al.(2001) 
presented a map of a line at an EME transition frequency 
towards SgrB2 of the order of 20". If we assume W51e2 to 
have this same source size the column density would have 
to be corrected to lower values by nearly an order of mag- 
nitude. More realistic, our best fit value (fit II) assumes 
an 5" size of W51e2. 



Furthermore, it can be seen from Fig. [7| (see simulated 
spectrum for N = lx 10 15 cm~ 2 ) that in some frequency 
regions strong EME lines, i.e. lines well above 3<r level, 
should reveal their internal line splitting. Due to the weak- 
ness of the observed lines no line was detected where the 
torsional splitting could be resolved so that there is no 
unique proof that these lines indeed belong to EME. 
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5.2. The other sources 

Sample spectra showing all examined sources at 245 
GHz are shown in Fig. |H1 At this frequency one of 
the strongest observed EME candidate lines appears in 
W51e2. However, no line above la appears in these spec- 
tra for any other examined source. Also for the other fre- 
quency bands we get similar results. For example, the well 
pronounced feature at 150845 MHz in W51e2 is only at 
the 1-2ct level in Orion KL and is completely invisible 
in G34.26 and NGC6334I. An upper limit for the col- 
umn density of G34.26, NGC6334I and Orion KL has 
been estimated by the use of myXCLASS and indepen- 
dently by RTDs using integrated noise at EME transi- 
tion frequencies and assuming an appropriate line width. 
Fig.0compares the 215 GHz region of W51e2 with Orion 
KL and also shows simulated spectra of EME for N= 
2 x 10 14 cm" 2 and N= 1 x 10 15 cm" 2 at 70 K (assum- 
ing an 5" source size). It is immediately clear that the 
column densities of these sources must be of the order of 
10 14 cm -2 or less rather than 10 15 cm" 2 as estimated by 



EME in W51e2. According to Table Hand the model by 



Charnley et al.(2001) 



Our results are summarized in Tabled The previously 
published line at the EME 53 2-42 3 transition frequency 



at 160.1 GHz in Orion KL by |Charnley et al. (2001)] has 
not been considered in the fit. 



6. Discussion 



Charnley et al.(2001) 



estimate the column density of EME in Orion KL and 
SgrB2(N) to be between 10 14 - 10 15 cm" 2 and in W51 
to be a few times 10 14 cm" 2 . EME column densities at 
around 10 14 cm" 2 correspond to signal strengths which 
are already close to the confusion limit in the HMC spec- 
tra. Nevertheless, long integration times and the availabil- 
ity of spectra at many frequencies enable us to estimate 
the column density of Orion KL to be below 1 x 10 14 cm -2 
and those of G34.26 and NGC6334I to be of the same or- 
der. 

We can now use our values in Table of 
W51e2 to draw some conclusions using the model of 
Charnley et al.(1995) According to Fig. 2 and Fig. 3a 
in their paper, which assumes a 10 -6 fractional abun- 
dance XeoI of ethanol as a starting value for the chemical 
evolution in their model, the age of W51e2 should be ei- 
ther > 10 5 yr if we use CH3OH or C2H5OH as a probe 
molecule or 10 4 - 10 5 if we look at CH3OCH3. However, 
that would imply a fractional abundance between < 10 -10 
and ~ 3 x 10~ 9 for EME. On the other hand, if we use 
Fig. 4a of ICharnley et al.(1995) with ethanol starting at 
XeoI=10~' this translates to Xeme < 10" 12 and X^me 
~ 4 x 10 -10 . Assuming C2H5OH to be the best molecule 
to determine the chemical age of the HMC it appears that 
in both cases, i.e. with ethanol starting values or Xe i 
= 10" 6 and XeoI = 10 -7 , the fractional abundance of 
EME is nearly the same, namely Xeme ~2x 10~ 10 . This 
value is close to our estimated abundance 5.6 x 10" 10 for 



Charnley et al.(1995)| Orion KL and NGC6334F should be 
roughly of the same age as W51e2. Although we have only 
upper limits on the EME abundance in these sources it is 
clear that as a consequence the C2H5OH fractional abun- 
dance at yr must be 10 -7 rather than 10 -6 . G34.26 is 
slightly older and SgrB2 is the oldest of the here discussed 
sources and should have an EME abundance Jeme of less 
than 10" 10 (corresponding to less than 10 14 cm" 2 column 
density as we see it for Orion KL, G34.26 and NGC6334I). 
This however is in contradiction to the claimed obser- 



vations of EME towards SgrB2 by |Charnley et al.(2001) 
when they estimate an abundance of 10 _iu - 10 _a . 

Looking at the observed values only, we notice that 
Orion KL has the highest fractional abundance of ethanol 
relative to H 2 followed by NGC6334I, G34.26, W51e2 and 
then SgrB2(N). However, no EME has been detected in 
Orion KL. The relative abundances of DME with respect 
to methanol and ethanol are shown in Fig.El Here, W51e2 
shows a relative high dimethyl ether content relative to 
methanol and ethanol compared to the other sources, 
which may be a general indication of a higher ether abun- 
dance relative to the alcohols. Assuming that we indeed 
have detected EME in W51c2 the ratio between EME 
and ethanol is 0.13 which is 5 times less than the ra- 
tio of DME and methanol. If all sources would have a 
0.13 EME/ethanol ratio the theoretical column densities 
would be around 2xl0 14 cm" 2 or higher 6 which at least 
for G34.26, NGC6334I and Orion KL has been shown to 
be not the case. Thus, W51e2 is indeed particularly rich 
in EME. 

We conclude that only the tentative detection of EME 



by Charnley et al.(2001) towards W51e2 can be confirmed 
within the following caveats 1) the EME Boltzmann fits 
include a few lines close to the noise level, 2) the obtained 
temperature of 70 K for EME is lower than that of other 
observed molecules (> 140 K), 3) no torsional splitting 
of the EME lines could be observed and also other char- 
acteristic spectral fingerprints have not been observed. 
However, there are lines at EME transition frequencies 
which until now can best be explained by the existence of 
EME in W51e2. 
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Source 




uec 


rLq. 




telescope 










(km s- 1 ) 




G34.26 


18 53 18.496 


01 14 58.66 


2000 


+60.0 


IRAM 30m 


NGC 6334 I 


17 17 32.300 


-35 44 04.00 


1950 


-8.0 


SEST 15m 


Orion KL 


05 35 14.171 


-05 22 23.10 


2000 


+8.8 


IRAM 30m 


W51 e2 


19 23 43.959 


14 30 34.55 


2000 


+55.0 


IRAM 30m 
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Table 2. Rotational temperatures, column densities and fractional abundances X relative to H2 of CH3OH, C2H5OH, 
CH 3 OCH 3 and C2H5OCH3 (EME) from our and other selected works. 



source molecules 




T rot [K] 


N t [cm" 2 ] 




Reference 


G34.26 






(1.8 ± 0.2) xlO 16 






CH3OH 


a-type 


368 ± 16 


4.5xl0~ 8 


MacDonald et al.(1996)| 






368 ± 16 


(1.6 ± 0.2) xlO 16 


4xl0~ 8 


Thompson et al.(1999)f 




b-type 


336 ± 14 


(1.9 ± 0.1) xlO 16 


4.8xl0~ 8 


MacDonald et al.(1996)| 




a+b type 


96 ± 17 


(2.6 ± 1.4) xlO 16 


6.5xl0~ 8 


Ikeda et al. (2002)1 




a-type 


140 ± 14 


2.2+ 4 - 7 xlO 16 


5.5xl0~ 8 


this work 


C2H5OH 




94 ± 27 


(3.5 ± 2.0) xlO 15 


8.8xl0~ 9 


MacDonald et al.(1996)| 






75 ± 37 


(1.7 ± 2.1) xlO 15 


4.3xl0 -9 


Ikeda et al.(2002)| 






152 ± 66 


6.0l 8 ; 4 xlO 14 


1.5xl0 -9 


this work 


CH3OCH3 




137 ± 32 


(3.7 ± 1.8) xlO 15 


9.3xl0~ 9 


Ikeda et al.(2002)| 






146 ± 38 


3.5+H xlO 15 


8.8xl0~ 9 


this work 


C 2 H 5 OCH 3 




100 - 150 


< 7 xlO 13 


<1.8xl0~ 10 


this work 


NGC 6334 I 




79±f 


3.1+S-g xlO 16 






CH3OH 




1.6xl0 -7 


Nummelin et al.( 1998b) | 






66 ± 11 


i.2±i:| x 10 16 


6.0xl0~ 8 


this work 


C 2 H 5 OH 




94+^ 


1.4+^3 xlO 15 


7.0xl0 -9 


Nummelin et al.( 1998b) | 






120 ± 44 


4.0+i' 2 x 10 14 
w — 0.9 


2.0xl0~ 9 


this work 


CEUOCHo 




20/440 


3.5xl0 15 / 1.3xl0 16 


4.1xl0~ 8 


Nummelin et al.( 1998b) | 






119 ± 39 


2.3tH x 1015 


1.2xl0 -8 


this work 


C2H5OCH3 




50 - 150 


< 8 x 10 13 


<4xl0~ 10 


this work 


Orion KL 






5 x 10 16 


5xl0~ 7 




CH3OH 




120 - 140 


Sutton et al.( 1985)| 






140 


4.5 xlO 16 


4.5xl0 -7 


Turner(1991)| (Onsala) 






192+| 8 

— ID 


(1.6±0.2) xlO 16 


1.6xl0 -7 


Turner(1991)|(NRAO) 






> 200 


2.5 - 30 xlO 16 


1.6xl0 -6 


Sutton et al.( 1995)1 






120-200 


4-8 x 10 16 


6xl0~ 7 


this work 


C2H5OH 




216 +3 f 90 

w — 105 

289 ± 330 


6.6+^ x 10 15 
4.8±l% x 10 14 


6.6xl0~ 8 
4.8xl0 -9 


Turner (1991) 
this work c 


CH3OCH3 




75 


2.5 x 10 15 


2.5xl0 -8 


Johansson et al.(1984)| 






63 


3.0 x 10 15 


3.0xl0 -8 


Blake et al.(1986)| 






9ll 3 i9 


1.3±o* xlO 15 


1.3xl0 -8 


Turner(1991)| 






125 


~ 3 xlO 15 


3.0xl0 -8 


Sutton et al.( 1995)| 






197.4l.tf 


l.l±g;g xlO 16 


l.lxlO -7 


Lee & Cho(2002)| 






126±43 


4.6±H x 1015 


4.6xl0 -8 


this work 


C2H5OCH3 




100-150 


<7 x 10 13 


<7xl0~ 10 


this work 


SgrB2(N) 






(2.8 ± 0.07) xlO 16 


5.6xl0 -9 




CH 3 OH d 


I 


49 ± 1 


Cummins et al.(1986) 






26±{f 


8. ll| xlO 15 


1.6xl0 -9 


Turner(1991)| 




II 


120 ± 20 


(3.5 ± 0.5) xlO 15 


7xl0 -10 


Cummins et al.(1986)| 






204 + o? 2 

w — 81 


6.9i3? xlO 15 


1.4xl0~ 9 


Turner(1991)| 


C2H5OH 




14 


5.3 xlO 14 


l.lxl0~ 10 


Cummins et al.(1986)| 






36 


4.8 + 2 4 xlO 15 


9.6xl0~ 10 


Turner(1991)| 


CH3OCH3 




24 


5 xlO 14 


lxlO -10 


Cummins et al.(1986)| 






801P 1 


2.3±° i xlO 15 


4.6xl0~ 10 


Turner(1991)| e 


W51e2 






1.0 xlO 17 






CH3OH, 




208 


2.7xl0 -7 


Ikeda et al.(2002)| 




a-type 


101±4 


2.5l°j xlO 16 


6.9xl0 -8 


this work 


C2H5OH 




169±28 


(3.1 ± 1.1) xlO 15 


8.6xl0" 9 


Ikeda et al.(2002)| 






143±20 


1.6+0- 4 xlO 15 


4.4xl0 -9 


this work 


CH3OCH3 




145±8 


1.6+° 4 xlO 16 


4.4xl0 -8 


this work 


C2H5OCH3 




69±94 


2+f 9 xlO 14 


5.6xl0 -10 


this work 



a We used N(H 2 )(G34.26) = 4xl0 23 cm" 2 , N(H 2 )(NGC6334I)= 2xl0 23 cm" 2 , N(H 2 )(Orion KL)= lxlO 23 cm" 2 , 

N(H 2 )(SgrB2(N))= 5xl0 24 cm" 2 and N(H 2 )(W51e2) = 3.6xl0 23 cm" 2 , see sectionH 

Observed at center of G34.26, as opposed to the halo. 
c The Onsala and CIT survey (see |Turner(199l"J| did not find any ethanol towards Orion KL and estimated an upper limit 

between 1.5 and 4xl0 14 cm -2 . 

|Cummins et al (1986)| 

Q and |Turner(1991)| used two rotational temperature diagrams (I and II) to fit the observed CH3OH lines. |Cummins et al.( 1986) 
uses the transitions with E u below 70 K for the fit I whereas |Turner ( 199 1 ) | splits the two components I and II at E u = 40 
K. 
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Transition 


Observed 
frequency 
(MHz) 


Theor. 
frequency 11 
(MHz) 


Obs-calc 
(MHz) 


line 
components 


(K) 


/ T mb dv 
(K km s" 1 ) 






llz.io-Hi.ii 


80 882.3(5)' 

SI 1 QK A(~\ \ 


80 882.7(10) d 

SI 1 QS *7(K\ 


-0.4 
— u.o 


5 




30.1 

1 1 S ^ 
1 lo.O 


0.09(4) 

U. ±0^0 ) 


□ 


u 

Cl 


353.32- 352,33 

363.33- 362,34 
042,32-0*1,33 
343,31-342,32 
1 3,34-o ( 2,35 


91439.6(5) 
y i oou. ± 

01 SI 9 7f A\ 


91440.4(10) 
91475.6(9) 
y 1 Dou.o^o ) 
91688.3(10) 

QI SI 9 

y 1 oiz.Q^y ) 


-0.8 

n 7 

n 1 
u. 1 


5 
5 

r. 
O 

5 

'J 


255.2 
269.3 

9^7 1 

241.5 

9S^ 7 
Zoo. / 


1.0(44) 
u.uo^y ) 


n 

1 1 

n 


b 
b 

C2, n 

11 

n 


352,33"35l,34 
293,26-292,27 
32,1" 2l,2 


Qf\ QQf) RCi\ 

yu oyu.o^o ) 
96464.1 6 


96 207.0(9) 
96464.3(6) 


u.u 

-0.2 


5 

5 


250.8 

1 7Q 1 

6.9 


0.002(5) 


□ 


b 

K 
1 ) 

11 


223,19-222,20 


107 654.8 6 


107 656.8(14) 


-2.0 


5 


108.3 


0.05(7) 


□ 


11 


72,5-61,6 

15l,15-14 ,14 


131350.2(10) 
131372.68(80) 


131350.7(8) 
131372.8(5) 


-0.5 
-0.12 


5 
5 


15.4 
46.7 


0.02(6) 
0.02(10) 


□ 
□ 


C3, n 
C4, n 


36l,35-36o,36 
344,30-343,31 
136,2-145,,, 
20o,20 - 19l,19 


150 661.3(6) 
150 794. V 


150 502.4(19) 
150 662.9(16) 
150 794.1(18) 
1 ^9^ 


-1.6 
0.0 
n 09 

U.Ui 


5 
5 
10 

K 
O 


260.06 
248.7 
76.6 

OU.'i 


0.12(2) 
0.17(0.05) 
n r >^\{^ 1 

u. £uy 1 1 j 


• 


b 
b 

C5 
Cfi 


[...]' 

iy5,x~^y4.i5 

[...]' 


91 ^ A e 


91 ^ ^!9fi ^fl ^ 

ill) OiiU.O^lO^ 


n q 


« 
U 


1 ri9 9 


u.^oyiu j 




b 
b 


28o,28"27l,27 


217 939.6(3) 


: 1 -| -7 /1 .1,1 7/1 \ 

217940.7(1) 


— 1.1 


5 


154.7 


3.66(19) 




b 


163,14-152,13 
31l,31"30o,30 


245 106.4(32) 
245 275.0(13) 


245 104.7(29) 
245 274.1(1) 


1.43 
0.86 


5 
5 


62.9 
188.8 


0.46(65) 
0.87(77) 


• 
• 


C7 


173,15-162,14 

282,27"27l,26 


252 190.7(11) 

253 309.3(3) 


252 189.4(29) 

253 308.2(12) 


1.3 
1.09 


5 
5 


69.5 
161.0 


0.70(79) 
16.01(126) 


• 


C8 

b 



a Rest frequencies (assuming a source LSR velocity of +55 km s _1 ) in W51e2. 

The numbers in parenthesis are the estimated uncertainties (la). 
c Frequencies are taken from |Fuchs et al.(2003)l 

d The deviation from the mean frequency of the transitions to its extreme values, e.g. due to rotational and internal rotational 
splitting the transition lla,io-lli,ii is split in 5 transitions which occur between 80 881.7 - 80 883.6 GHz, thus the frequency 
denoted in the table is 80 882.7(10). 

e Manually determined (or in case of noise the closest peak). Here the transition 136,x-15s, y are with k c = 7-9, 7-10, 8-9, 8-10 
and the x in 195,2,-194,15 denotes either 14 or 15. 

f Region with many blends of EME with other (dominating) molecules. 

• Lines included in the Boltzmann fit II, see Fig. 
□ Lines included in the Boltzmann fit I and II, see Fig. |S] 

n "Noisy" line, i.e. a line which has either less than 2a intensity or where the intensity is difficult to obtain due to baseline 
variations. 

b Denotes a blend with an other molecular line where the other line contributes most to the total intensity. 
Ci Denotes a contribution of another molecule to the total intensity where the blended line is estimated to contribute less 
than 20% (if not otherwise mentioned) compared to EME. Assuming T ro t ~ 100 K: 

ci C 2 H 5 OOCH with TV <5-10 15 cm" 2 (could be up to 50%), c 2 HOONO2 with TV <5-10 14 cm" 2 , c 3 C-H2C3O with TV <2-10 12 
cm" 2 , c 4 C2H3NH2 with TV <8-10 15 cm" 2 , c 5 C3H7CN with TV <2-10 14 cm" 2 , c 6 C 2 H 5 C 15 N with TV <3-10 13 cm" 2 based 
on our observation that C2H5CN has TV ~ 3-10 14 cm" 2 , c 7 CH 2 (OH)CHO with TV <8-10 14 cm" 2 , c 8 a-H 2 C=CHOH with 
TV <2-10 14 cm" 2 . 
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Table 4. Estimated column density for various source sizes and rotational temperatures for W51e2. 



source size 


Trot 


column density 


[arc sec] 


[K] 


[cm- 2 ] 


20 




c , i a13 

~ 5 x 111 




ou 






70 


- 4 x 10 13 


10" 


50 


- 8 x 10 13 




70 


- 9 x 10 13 




150 


1.3 x 10 14 


5" 


70 


2.0 x 10 14 




150 


4.5 x 10 14 
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Table 5. Observational parameters 
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Receiver 
frequency 
(MHz) 


rp(a) 

-tsys 

(K) 


T 


Tot. integr. 
time (on+off) 
(min) 


p(b) 


Telescope 
beam HPBW 

(") 


lcr (rms) 
noise level 
(K as T£) 


Band 
width 
(MHz) 


Frequency 
resolution 
(MHz) 


G34.26 


















81 041 


144 


0.03 


85 


0.79 


30 


0.012 


490 


0.3125 


91 626 


122 


0.03 


115 


0.77 


26 


0.010 


490 


0.3125 


C\0 OT1 

96 371 


1 O T 

137 


A AO 

0.02 


1 /ir 

145 


A TT 

0.77 


A C 

25 


A Al A 

0.014 


A AA 

490 


A OI OC 

0.3125 


1 oi OTO 

131 372 


200 


A A/" 1 

0.06 


115 


A T1 

0.71 


18 


A Al O 

0.018 


A AA 

490 


A OI OC 

0.3125 


150 750 


OAT 

307 


A A A 

0.04 


220 


A £ O 

0.58 


15 


A Al T 

0.017 


A AA 

490 


A OI OC 

0.3125 


Oi r ran 

215 500 


088 


A OO 

0.22 


195 


a c a 

0.5b 


11 


A AOO 

0.0o2 


1 AAA 

1020 


1 A 

1.0 


217 900 


OCA 

o50 


A 1 A 

0.14 


1 A T 

14/ 


A C C 

0.55 


11 


A AOO 

0.02,3 


1 AAA 

1020 


1 A 
1.0 


245 260 


r 1 1 

511 


1 A A 

1.90 


1 ^n 

160 


A C A 

0.50 


1 A 

10 


A AO A 

0.034 


CT 1 A 

510 


1.25 


252 227 


O A O 

848 


A O T 

0.37 


85 


A AO 

0.48 


1 A 

10 


A A ET A 

0.059 


CT 1 A 

510 


1.25 


lvr /" 1 />•)•> < i 

JNCjtC 6334 1 


















o 1 1 no 

81 198 


OAT 

327 


0.41 


O O 

38 


A TC 

0.76 


CT A 

59 


A AAA 

0.009 


AA A 

994 


a c m 1 

0.6911 


91 450 


195 


0.26 


37 


0.74 


55 


0.007 


994 


0.6911 


A 1 /?OA T 

91 630.7 


185 




87 


0.74 


55 


0.005 


994 


0.6911 


91 811.7 


A1 A 

212 


A A/" 1 

0.06 


r T 

57 


A T A 

0.74 


r r 

55 


A AAT 

0.007 


AA A 

994 


a c m 1 

0.6911 


a/"* oon 

96 380 


252 


A A 1 

0.21 


55 


A TO 

0.73 


C O 

53 


A AAA 

0.009 


AA A 

994 


A c m A 

0.0914 


150 662 


A/" 1 o 

263 


A O A 

0.34 


A CT 

95 


A £ CT 

0.65 


O O 

33 


A Al A 

0.012 


AA A 

994 


A £Z AD A 

0.6989 


A /i CT 1 AC 

245 11)5 


513 


A A/" 1 

0.96 


TT 

77 


A A T 

0.47 


22 


0.024 


AA A 

994 


A £Z AO A 

0.6989 


A /i P* A T /) 

245 274 


534 


0.50 


127 


0.47 


22 


0.017 


994 


A /"* A O A 

0.6989 


Urion KL 


















O 1 A /f 1 


1 O O 

138 


A AO 

0.03 


i oc 

125 


A TA 

0.79 


OA 

30 


A A 1 A 

0.012 


A AA 

490 


A OI OC 

0.3125 


91 626 


1 O A 

182 


A AT 

0.07 


A c: 

95 


A TT 

0.77 


26 


A A 1 A 

0.014 


A AA 

490 


A OI OC 

0.3125 


i rn tc a 

150 750 


O 1 o 

316 


A 1 A 

0.14 


AAA 

220 


A C O 

0.68 


16 


A AAA 

0.029 


A AA 

490 


A OI OC 

0.3125 


215 500 


CT O A 

539 


A O O 

0.33 


1 ^n 

160 


A CT 

0.56 


-i -i 
11 


A ATO 

0.078 


1 AAA 

1020 


1 A 

1.0 


245 275 


740 


0.39 


190 


0.50 


10 


0.053 


510 


1.25 


W51e2 


















O 1 a A 1 

81 1)41 


1 1 O 

113 


A A A 

0.04 


A c: 

95 


A TA 

0.79 


OA 

30 


A A 1 A 

0.012 


A AA 

490 


A OI OC 

0.3125 


91 525 


lip 

lib 




60 


A TT 

0.77 


A/^ 

26 


A Al A 

0.012 


A AA 

490 


A OI 

0.3125 
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Fig. 1. Stick spectrum of EME at a rotational temperature of 130 K which is typical for complex molecules in 
hot core regions. EME has a dense spectrum with many characteristic features. The frequencies are taken from 
Fuchs et al.(2003) The (*) indicate our chosen frequency regions for the W51e2 obseravtions. The intense line at 237 
GHz was not selected because of a blending line at the same frequency. 
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Fig. 2. Observed spectrum of W51e2 at 151 GHz (upper) and simulated spectrum using myXCLASS (lower) including 
C2H5OCH3, C 2 H 5 OOCH, CH3CH3CO, a-(CH 2 OH) 2 , c-C 2 H 4 0, C 2 H 5 OH, CH3OCH3, CH 3 CHO-e, CH 3 OH, C 2 H 5 OH, 
CH 3 OCHO-e, and CH 3 OCHO-a . 
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Fig. 3. Spectra recorded towards W51e2. The EME transitions 20 ,20-19i,i9 at 150 845.5 MHz, 19 5iX -19 4i i 5 at 215 325.4 
MHz , 16 3i i4-15 2 ,i3 at 245106.4 MHz , and 17 3 ,i 5 -16 2 ,i4 at 245190.7 MHz are indicated. 
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Fig. 4. The 20 ,20-19i,i 9 transition of EME at 150 845.5 MHz. The line profile observed towards W51e2 (bold) is 
compared to three simulated line profiles (dashed) of 0.2, 0.7, and 1.4 MHz line widths (FWHM). For a line width 
of 0.2 MHz (FWHM) the splitting due to internal torsional motion of the molecule is partly resolved. The stick 
diagram indicates the position and relative intensities of the split lines based on our recent high resolution laboratory 
measurements Fuchs et al.(2003) For comparison with the observed spectrum all simulated spectra are normalized to 
have the same peak intensity. 
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Fig. 5. Boltzmann plot for the EME transitions found towards W51e2 based on 14 observed lines, assuming no beam 
dilution. Due to difficulties in assigning the lines to EME two fits have been been performed. I) Fit including all lines 
marked with □ and • in TableEl Here T r ot — 126±115 K and log(Nt)=14.9±0.7. II) Fit using only lines with intensities 
> 2<T (•) resulting in T ro t — 54±38 K and log(Nt)=14.3±1.2. Lines which are less than 2a intense are marked with 
an □. 
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Fig. 6. Observed spectrum at W51e2 (top, dashed) and calculated spectra for different column densities of EME at 
70K. 




215300 215400 215500 215600 

frequency [MHz] 



Fig. 7. Observed spectra of Orion KL and W51e2 at 215 GHz. This frequency region has the highest line density of 
all our observed bands but even here gaps can be found where the presence of EME can be checked. Two simulated 
spectra of EME are plotted, i.e. N=2 x 10 14 cm" 2 (blue) and N=l x 10 15 cm" 2 (red) with T rot =70 K. The splitting 
of the lines which can best be seen above 215600 MHz is mainly due to the internal rotation of the methyl groups 
of EME with an additional small asymmetry contribution, a) For Orion KL the previously estimated limits of 10 14 
- 10 15 cm -2 by Charnley et al.(2001) seem questionable as can be seen from the superimposed simulated spectra on 
the observed spectra, b) In the case of W51e2 the simulated spectrum is in agreement with the observation. 
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Fig. 8. Spectra of the observed sources at 245 GHz. Only W51e2 reveals a line at an EME transition frequency, 
whereas all other sources do not show any significant feature above la level. However, W51e2, G34.26 and NGC6334I 
have similar spectra indicating the same chemical composition and equivalent dynamics. 
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Fig. 9. Relat ive abundance s of CH3OH, C2H5OH, CH 3OCH3 derived from Table |3 by using the values from 
Turner(1991)| for Orion KL, |Turner( 1991 )| for SgrB2(N), |Thompson et al.(1999)| and |Ikeda et al. (2002)1 for G34.26, 



this work for W51e2 and |Nummelin et al.(199 8b) for NGC6334I. The ratio of EME/C 2 H 5 OH in W51e2 is 0.13 which 
is 5 times smaller than the ratio of CH3OCH3/CH3OH. 



